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Using resonant interaction between atoms and the field in a high quality cavity, we show how to 
realize quantum random walks as proposed by Aharonov et al [Phys. Rev. A 48, 1687 (1993)]. The 
atoms are driven strongly by a classical field. Under conditions of strong driving we could realize an 
effective interaction of the form iS x (a — a') in terms of the spin operator associated with the two 
level atom and the field operators. This effective interaction generates displacement in the field's 
wavefunction depending on the state of the two level atom. Measurements of the state of the two 
level atom would then generate effective state of the field. Using a homodyne technique, the state 

in : of the quantum random walker can be monitored. 

O 

o 

In a very interesting paper Aharonov et al Q] proposed the idea of a quantum random walk. Here a random walker 
is constrained to move left or right depending on the state of an auxiliary quantum mechanical system. One then 
examines the state of the random walker subject to the measurement of the state of the auxiliary system. As an 
interesting consequence of this quantum random walk, Aharonov et al 1] found that the walker's distribution could 
shift by an amount which could be larger than the width of the initial distribution. Further the displacement could 
, be much larger than the classical displacement. Several proposals QHIESIEIESSl exist for realizations of the 
quantum random walk. For example Aharonov et al gave a cavity QED model where the photon number distribution 
can get displaced. Sanders et al Q considered a dispersive interaction in the cavity of the form S z (a + a') and 
considered the random walk of the field on states on a circle. Other interesting theoretical schemes for implementing 
quantum walks have been suggested in ion-traps 3 and in optical lattices [j| . Knight et al @ further showed that an 
earlier experiment Q was a realization of quantum random walks. A scheme using linear optical elements has been 
' recently implemented §]. 

Here we propose a method which yields precisely quantum random walk as proposed by Aharonov ct al. We use 
-4— > ' cavity QED however we drive the atoms by an external field. Currently there is considerable progress in realizing a 
variety of high quality cavities and a variety of interactions and thus one is in a situation where proposals like the 
one presented here are likely to be implemented. 
q-( The organization of the paper is as follows. In Sec. II we present the details of our model and show the conditions 
• • ■ under which such a model gives rise to an effective Hamiltonian which we use in Sec. Ill to realize quantum random 
walk. In this section we also present the results for the Wigner function for the state of the quantum walker. In 
Sec. IV we show how the homodyne measurements of the field can be used to check the characteristics of the quantum 
random walk. In Sec.V we incorporate the effects of decoherence due to the decay of the field in the cavity. In the 
appendix we discuss the state of the walker if no conditional measurements are made and establish relation to classical 
random walks. 



II. EFFECTIVE HAMILTONIAN FOR QUANTUM RANDOM WALK USING DRIVEN ATOMS 

We consider a two level Rydberg atom having its higher energy state |e) and lower energy state \g), interacting with 
a single mode of the electromagnetic field in a cavity. The atom passes through the cavity and interacts resonantly 
with the field. Further the atom is driven by a strong classical field. For simplicity we choose atomic transition 
frequency the cavity frequency and the frequency of the driving field to be same. The Hamiltonian for the system in 
the interaction picture is written as 

H = -ihg (S+a - a* S~) + h (S + £ + S~S*) , (I) 
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where g and £ are the coupling constants of the interaction of the atom with the cavity field and with the deriving 
field. We have chosen g as real and £ as complex. The annihilation (creation) operator for the field in the cavity is 
a(<r) and S + , S~ are atomic spin operators. The last term in Eq.QJ is the interaction with the external field. We 
further rewrite the above Hamiltonian in a picture in which the interaction with the external field has already been 
diagonalized. 

\i>) =e lht \iP); h = S+£ + S-£\ (2) 
where \ip) is transformed atomic state in new picture from old atomic state \ip). The Hamiltonian in this picture is 



H = -ige lht {S+a - S-^)e~ lh \ (3) 
ih 

W\ 



e iM = cos(\£\t) + ^sm(\£\t). (4) 



The atomic spin operators S transform as 



\£\ 



■^S* S m(\£\t)co S (\£\t), (5) 
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e 



iht S - e -iht = g- cos 2Qg| t ) + £ sin 2(|£|t)5+ 



2 -^rS z sin{\£\t)cos{\£\t). (6) 



Using Eqs.Q and ©, Eq.Q becomes 



H = -ig (Vcos 2 (|£|i) + ^ s in 2 (\£\t)S- 

04 f* \ 

— \£\ SZ sin ^ £ C0S ^ '7 a ~ H ' c ' ^ 

We note that the Hamiltonians of the above form have been previously used to treat the inhibition of the spontaneous 
emission [j| and for the production of Schrodinger cat states [ljj . We assume that the atom is driven strongly so that 
\£\ is large and hence we drop rapidly oscillating terms from Eq.Q i.e. e ±2l l £ l* =>. 0. Then Eq.(J7J) reduces to 

H = - % 4 (s+ + ^S-)a-H.c. (8) 



\£\ 

We choose £ * 2 /\£ | 2 = 1, in general, this can also be done by adjusting phases with atomic operators. Then the Eq.© 
takes the form 



H eff =gS*{-—y (9) 

Note the appearance of the well known displacement D(a) = (a^a — aa*) in the Eq.Q. In particular we have the 
momentum operator (out of phase quadrature for the field). Further it should also be noted that h as defined by 
Eq. commutes with H e f f . In the original interaction picture the Hamiltonian for our model will be 

H eff =gS*(±^-)+2\S\S*. (10) 



In the effective Hamiltonian l|l(J|l field displacement operator appears with atomic operator, which can produce dis- 
placement in field state depending on the atomic state. 
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III. REALIZATION OF RANDOM WALK 



We next examine the evolution of the system of the two level atom and the field inside the cavity. Let us consider 
that, initially the atom is in the superposition state |<&) = (ci|e) +c 2 \g)) and the field is in a coherent state \a). Using 
Eq. HI U| l the combined state of the atom-cavity system after time t is given by 



W)) 



2i\E\tS x ] |$)|a) 



exp [gtS x ((J - c 

-(\g) + \e))\a + gt/2) 



c+e_ 
2 



(\9) 



\g) 



+ 



\e))\a-gt/2), 
3 + gt/2) 
* + gt/2) 



2 



-\a-gt/2) 
■\a-gt/2) 



= (\£\ + g -Im{a)) t; 



(11) 
(12) 

(13) 
(14) 



where c+ = c\ + c 2 and c_ = c\ — c 2 . Using normalization of atomic states we can select c_/c + = tan#. Thus the 
detection of the atom in state |e) or \g) leaves the cavity field in a superposition of states \a + gt/2) and \a — gt/2). 
For small values of gt the states \a + gt/2) and \a — gt/2) overlap completely and thus quantum interference effects 
between \a + gt/2) and \a — gt/2) becomes significant. If we assume that the atom is detected in its ground state \g). 
Then the state of the field inside the cavity can be written as 



\ijj f ) tx \e' mt D{gt/2) + e l|£|t tan(6»)L>(- 5 t/2) 



(15) 



Clearly after passing one atom through the cavity the field inside the cavity is displaced backward or forward along 
the line in a random way by the step of gt/2. We can now iterate the above step to obtain the state of the field after 
the passage of N atoms. We assume that atoms enter in the cavity in the state |$) and after interaction with the 
field inside the cavity detected in their ground state \g). Note that the displacement operators appearing in the above 
state commute each other [D(gt/2), D(— gt/2)] = for real gt. Thus the field state after the passage of N atoms is 
given by 

N 



\^ f (N)) = C \e- mt D(gt/2) + e mt t&n{9)D{~gt/2)\ [ \a) 



N 
m=0 



im\£\t jjm 



J(N- m )\£\t 



(tan( 



\N — mjjN—m 



N 



= C J2 

D 



m=0 
N—2m 



i(N-2m)\£\t 



(tan 9) 



.a 

N-r 



a 



(-fl*/2)|«), 



N 



J(N-2m)4> 



(tanfl) 



N—m 



m=0 

\a-(N-2m)gt/2)), 



(16) 



(17) 



where C is normalization constant and we have used the property of the displacement operator D (ex) — D(—a). 
On writing the above result in coordinate space representation, we get the wavefunction ipN(x, ct) = (x\ipf(N)) 



4>n(x, a) = C 



n 

E 

m=0 



i(N-2m)(p 



(tan 9) 



N- 



ip a (x + [N-2m]l), 



(18) 
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FIG. 1: The probability distribution P(x) for the position of the quantum random walker, assuming initial wave packet as 
Gaussian exp[— (x — a) 2 /2] for a = 0, step size I = 0.05, (j> — 2ty and 6 — 2tv/3. 

where ^p a {x) = {x\a) is the wavefunction corresponding to the initial cavity field state \a) which is centered at x = a 
and the step size of the random walker is I = gt/2. We note that we have recovered the result of Aharonov et al 
In Fig^we have plotted the probability amplitude distribution for initial wave function ip a { x ) ~ exp[— (a: — a) 2 /2] 
for real values of x and a = 0. The displacement depends on 9, <p and the number of steps N . The unexpected 
displacement in the state of the random walker is the result of constructive quantum interference between the states 
generated in various steps which comes from the off diagonal terms in P(x) = |^>jv(a, x)\ 2 . We have checked this by 
dropping the off diagonal terms in P(x), in that case P(x) remains same in shape as the initial wave packet but shifts 
by an amount Nl. The displacement of the random walker is not bounded by the classically possible maximum and 
minimum displacements ±Nl. The quantum interference leads to an arbitrary displacement in the random walker's 
position and can be much larger than ±Nl. A small squeezing in wavepacket is also generated from these interference 
effects. The selection of phase <fi is also critical for displacement in the position of quantum walker, for example for 
the parameters used in Fig^ the maximum displacement in the position of quantum walker occurs when <f> is integer 
multiple of 7r and there will be minimum displacement when <fi is half integer multiple of n. 

For visualizing quantum interferences we plot the Wigner function of the random walker in Fig|2] The Wigner 
distribution for any state ip(x) can be obtained by using the definition 



In the Fig|2Ia) the field is in its initial coherent state and the wigner function is perfect Gaussian. As the field is 
displaced by random steps, by passing atoms through the cavity, quantum interference effects start deforming the 
shape of the Wigner function from the Gaussian. After few steps the Wigner function is squeezed in x quadrature 
and gets displaced by an arbitrary distance in x. In Fig^Jb), (see also FigQJa)), we have shown the Wigner function 
after 10 random steps for initial Gaussian wave packet. The squeezing is also clear from the Fig^ which shows the 
narrowing of the distribution P(x). It is clear that the displacement in the position of random walker comes as a 
result of quantum interference which is consequence of quantum coherence between the states generated in random 
steps. 



We next discuss how we can probe the quantum state of the random walker. We propose homodyne techniques |12( 
for measuring the state of the random walker. Such homodyne measurement can be performed by mixing an external 
resonant coherent field to the cavity and then probing the resultant cavity field by passing a test atom through the 
cavity. In the previous section, we have shown how the cavity field is displaced backward or forward in a random step 




(19) 



IV. MEASUREMENT OF THE STATE OF THE RANDOM WALKER 
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FIG. 2: The Wigner function W(x,p) of the state of the random walker, after number of steps (a) N — (b) N = 10, using 
same parameters as in FIG0 

by passing single atom through the cavity. The state of the field in the cavity after such N steps can be monitored 
by homodyne measurements which can be implemented in the same experimental set up. After displacing the field 
inside the cavity by N random steps, by passing N atoms, a resonant external coherent field \j3) is injected into the 
cavity. After adding the external field, the state of the resultant field in the cavity is 




m— 



D(fl)\a 



(N-2m)gt/2)), 




n 771 — 



(n\D({3)\a - (N - 2m)gt/2))\n), 
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FIG. 3: The probability of detecting probe atom in its ground state as a function of 5 for the state of the quantum random 
walker after number of steps N = ( solid line) and N = 10 (dashed line). The parameters used are same as in Fig^and the 
interaction time for the probe atom is selected such that gt p = 1.5tt. 



N , s 
m=0 ^ ' 



N-m 



(n\D(f3)\a-(N-2m)gt/2)). 



(21) 



Now we bring a similar atom in its lower energy state \g) to probe the cavity field. The probability of detecting the 
probe atom in its lower state \g) after crossing the cavity in time t p is 



The interaction time t v for the probe atom is selected such that if there are photons in the cavity it leaves the cavity 
in its higher energy state |e) with larger probability. If we choose the external field \(3) such that (3 = —a + 8, the 
probe atom will leave the cavity in its ground state with larger probability when the value of 6 will be opposite and 
equal to the displacement of the random walker from the initial position a. Thus the probability of the probe atom 
leaving the cavity in its lower state \g) would, as a function of <5, have peak corresponding to the positions of the 
random walker after TV steps. In FigEl we plot the probability of detecting the probe atom in its lower state with 
8. The solid line curve is result of homodync measurement of the position of the random walker corresponding to 
its initial state. The dashed line curve is corresponding to the homodyne measurement after 10 steps using the same 
parameter as in FigQ] Clearly the homodyne measurement yields the state of the quantum walker (Fig^). Thus the 
homodyne measurement can be an elegant way for monitoring the position of the random walker in our model of 
realizing quantum random walks. 



Quantum random walks are different from the classical random walks in the sense of quantum interferences which 
may lead much larger displacements in the position of quantum random walker than the classically possible maximum 
displacements. These quantum interferences are the consequences of coherence in the system. Clearly we need the 
coherence to live for a long time and thus it is important to study the effects of the decoherence of the system. In 
this section we study the decoherence of the state of the random walker due to damping in the cavity. This can be 
done using the master equation 




(22) 



n 



V. DECOHERENCE OF THE GENERATED STATE OF THE RANDOM WALKER 



P = 



— (a^ap — 2apa ! + pa' a) 



(23) 



where k is cavity field decay parameter and we carry analysis in the absence of thermal photons. For initial state Ijl7|) 
we find the density matrix after time t 
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FIG. 4: The decoherence of the state of the random walker in terms of Wigner function at different times, (a) for Kt = 0, (b) 
for Kt = l/4N 2 l 2 , (c) for Kt = l/2N 2 l 2 , (d) for Kt = 2/N 2 l 2 , other parameters are same as FIGH(b). 



N N , s , 
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m— n— 

(a - (JV - 2m)2|a - (N - 2n)/)) (1 ~ e_Kt) 

|a - (JV - 2m)0t(a - (N - 2n)l\ t , (24) 



where \() t = |Ce _Kt ' 2 ). In the limit Kt << 1 the Eq. simplifies to 



N N 



m—0 n—0 



Pit) = ICfY.YA 1 ) ( " )e^-^{Un6f N - 



e - 2Ktl2( - n ~ m)2 \a — (N — 2m)l)(a - (N - 2n)l\. (25) 

Thus the coherence of the state decays on the time scales \/2N 2 l 2 . In Fig0]we show the decoherence effects due to 
the cavity damping in the state of the quantum random walker in terms of Wigner function. As the time progresses 
from (a) to (d) the decoherence reduces the quantum interference effects and the state of the random walker decays to 
its initial state. In Fig^Ja) the Wigner function for the state of the random walker after 10 steps using the parameters 
of Fig[2Ib) is plotted which is squeezed in x quadrature and centered at x w —2. As a result of decoherence due to 
cavity damping the quantum interferences start decaying and the Wigner function changes to the perfect Gaussian 
shape, FigQJc) centered at x = Nl. Now the field inside the cavity is almost in coherent state and decays with the 
cavity damping rate. Further the life time for the state of the quantum random walker is given by Tpj = T c /2N 2 l 2 
where T c = 1/k is life time for field in the cavity. 
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VI. CONCLUSIONS 



In conclusion we have shown a simple possible realization of quantum random walks using cavity QED. We have 
proposed homodyne detection for monitoring the position of the random walker. We have also discussed the deco- 
herence effects and the time scales at which quantum nature of random walks survives. As a result of new emerging 
technologies various improved cavities are feasible these days which makes our proposal much interesting and 
realistic. Such realization of quantum random walks may be useful for implementing various algorithms [T^ based 
on quantum random walks. Finally it should be noted that the generalizations of the present work to more than one 
dimensions are possible. 



APPENDIX: STATE OF THE WALKER FOR NO MEASUREMENT ON THE ATOMIC STATE 

In this appendix we would like to connect the result Q13[) explicitly to the case of classical random walk. For this 
purpose we find the reduced state of the field from (|13|l . We also set Ci = 1, C2 = 0, then the reduced state of the 
field pf is 

P/ = 5(|a+f><a + fl + |a-f><a-f|) (A.l) 
Clearly the state of the field after the passage of N atoms would be 




\a + !L.(N-2m))(a + ^(N-2m)\ (A.2) 

fl\ N ±^ N\ gt ., gt , , A s 

= U E ES |a + i p)(ft+ i 1 ' 1 (A - 3) 

^ ' p=-N 2 - 2 

which is reminiscent of the result for classical random walk in the sense that the weight factor of the state \a + 
YP)(a + i s same as the probability of finding the walker at the site p It should however be borne in mind 
that the coherent states \a + ^p) and \a + ^p') are not orthogonal for p ^ p' . 
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